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We construct a novel observable for finite temperature QCD that relates confinement and chi- 
ral symmetry. It uses phases as boundary conditions for the fermions. We discuss numerical 
and analytical aspects of this observable, like its spectral behavior below and above the critical 
temperature, as well as the connection to chiral condensate, center symmetry and the canonical 
ensemble. 
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1. Motivation and introduction 

The main phenomena in QCD at the finite temperature transition are deconfinement and chiral 
symmetry restoration. There has always been the question whether there is a mechanism connec- 
ting the two, in particular since in the quenched theory the corresponding phase transitions occur 
at the same critical temperatures [1]. Here we define an observable that indeed links the two and 
discuss numerical findings from quenched lattice configurations [2] . 

The Polyakov loop (with denoting path ordering) 

P{x) = ^exp (i l^^dxoAo{xo,x)'^ , jS = l/kBT, (1.1) 

is the order parameter of confinement, being traceless in the confined phase and moving towards 
the center of the gauge group, for SU(3) the three elements {l,exp(2;r//3),exp(47r//3)} I3, at 
temperatures above Tc. This behavior can be understood by its relation to (exp(— jSf ) of) the free 
energy f of a single heavy quark, which is infinite in the confining regime. 

The spectral density of the Dirac operator at the origin, p (0) , on the other hand, is the order pa- 
rameter of chiral symmetry. It is related to the condensate by the famous Banks-Casher relation [3] 
(YY) = — ^P(O) and vanishes above Tc. 

How does confinement leave a trace in the Dirac spectrum? After all, the quarks should not 
only know about chiral symmetry, but also that they are (de)confined. The answer will lie in the 
dependence on temporal boundary conditions, as we will show now. 

2. Idea, derivation and interpretation of the new observable 

We use the lattice as a regulator. The (untraced) Polyakov loop 

No 

P{x) = Y[Uo{^,x), (2.1) 
is built from temporal links. For the lattice Dirac operator we use the staggered one [4] 

^(^^y) ^ ^ I^mW [U^{x)d,+^^y-h.c.] , r^^ix) = , (2.2) 

which can be viewed as hopping by one link. 

It is obvious and well-known that the ^-th power of the Dirac operator at the same argument, 
D''{x,x), contains all products of links along closed loops of length k, starting and ending at x. The 
Polyakov loop is among these loops for k = Nq, but how to distinguish it from 'trivially closed' 
loops (like, e.g., the plaquette), that do not wind around the temporal direction? 

The tool for this has been introduced by one of us in [5]. One needs phase boundary conditions 
for the fermions 

Y{xo + P,x) = e''PY{xo,x) , (p G [0,27r] . (2.3) 

The physical case of antiperiodic fermions is obtained for (p = n. These boundary conditions^ 
amount to an imaginary chemical potential. They can be easily implemented by replacing the 

'Note that fermion bilinears like xjf^ iff are strictly periodic. 
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temporal links Uq in some time slice by e"''Uo and likewise Uq by e^'fU^. As a consequence, all 
Polyakov loops get a factor e"f, inverse Polyakov loops a factor e^"^ and those with higher winding 
number get the corresponding power of that phase factor, while the trivial loops stay the same. 

In this way, the Polyakov loop can be reconstructed from the Dirac spectrum by using at least 
three boundary conditions, see [6]. This 'thin' Polyakov loop, however, has poor renormalization 
and scaling properties and it turned out that in this approach it is UV dominated [6]. 

Influenced by the Jena group [7], we instead consider the propagator with some probe mass m. 
The Dirac operator at a particular boundary condition (p is denoted by D(p and we use a geometric 
series to represent the propagator, 

This representation obviously contains all powers of the Dirac operator. Plugging in the definition 
(2.2) and the factors of e"'', the propagator is given as a product of links along all closed loops, 

where |/| is the length of the loop and sign(/) comes from the staggered factor. The ordered product 
of U^{x) is over all links (jc,/i) in the loop. 

Of importance in (2.5) is the phase factor, where counts how many times the loop winds 
around the temporal direction. One can project onto a particular winding ^ by a Fourier transform 
w.r.t. (p, 

1 r^^ 

d(pe-''P'>. (2.6) 



2k Jo 

Specifying to a single winding, q = I, like for the Polyakov loop, we arrive at [2] 

"^--'■^-(tr^— ) = -^ E ^(tr.n W). (2.7) 

271 V\ m + D^/ mV^(^j(2am)l'l \ ( ' / 

This completes the derivation of our new observable £ , which we refer to as the 'dual condensate', 
because it is obtained through a Fourier transform from the trace of the propagator. Indeed, in the 
massless limit (after the infinite volume limit as usual) we obtain the chiral condensate 

lim lim £ = - f^"" ^ e-"P lim lim (W)<p= — ^'""^P(0)(p , (2.8) 

m^ov^=» Jo 2k m^OV^^ ^ Jo 2 '-^ 

integrated with a phase factor over the boundary conditions. Making use of the Banks-Casher rela- 
tion at every individual angle (p, we furthermore obtain the representation in terms of the eigenvalue 
density p(0),p. 

The right hand side of (2.7) represents the 'dressed Polyakov loop', that is the set of all loops 
which wind once around the temporal direction. In the infinite mass limit, detours become sup- 
pressed and only the thin, straight Polyakov loop survives as it is the shortest possible loop in this 
set. 

We would like to stress that Eq. (2.7) is an exact relation and is valid for individual configura- 
tions. 
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Figure 1: The dressed Polyakov as an order parameter. Left: Expectation values for m = 100 MeV as a 
function of T at various lattices, i.e., for different volumes and lattice spacings. Right: The corresponding 
plot for the conventional Polyakov loop. 

3. Numerical results and more interpretation 

In the following we discuss various aspects of the relation (2.7). First of all, Fig. 1 shows that Z 
is indeed an order parameter. Keeping the mass m fixed, £ vanishes below the critical temperature 
(which is about 280 MeV in the quenched case) and develops an expectation value for higher 
temperatures. One finds that the results (when expressed in physical units) obtained for different 
volumes and with different resolution essentially fall on a universal curve. This illustrates the good 
renormalization properties of our observable, which are inherited from the renormalization of the 
chiral condensate. The improved renormalization properties may also be understood as an effect 
of the dressing which renders the new observable less UV dominated. For comparison we have 
plotted the corresponding expectation values of the conventional thin Polyakov loop in the right 
hand side panel of Fig. 1 . 

In a spectral representation, the Dirac operator D^p in Z can simply be replaced by a sum over 
all its eigenvalues, 

where the index (p on the eigenvalues Xjp^ again refers to the boundary condition angle. This is 
actually a suitable representation to numerically compute the dressed Polyakov loop (in contrast, 
the right hand side of (2.7) contains an infinite sum over loops even on a finite lattice). On our 
lattices we have calculated all eigenvalues and approximated the (^-integral by the trapezoidal rule 
with 16 equidistant boundary conditions. 

As the eigenvalues appear in the denominator, we expect the sum to be dominated by the IR 
modes. As Fig. 2 shows, this is confirmed by the lattice data, if m is not too large. 

How is a finite resp. vanishing order parameter £ built up by the eigenvalues? Fig. 3 shows 
that they respond differently to the boundary conditions in the confined vs. deconfined phase. In 
that figure we plot the expectation value of the propagator, i.e., the integrand of £ without the 
Fourier factor. The eigenvalues are independent of the boundary condition in the confined phase. 
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Figure 2: Individual and accumulated contributions to the spectral sum (3.1) at two different masses. 

which leads to a vanishing order parameter Z. In the deconfined phase, on the other hand, the 
eigenvalues show a typical cosine-type of modulation. Together with the Fourier factor this yields 
a nonvanishing £ (proportional to the amplitude of the modulation). 

The chiral condensate has to behave essentially in the same way, as it is the integrand in the 
massless limit. Although the chiral condensate is finite in the confined phase, it is independent of 
the boundary condition <p and hence results in a vanishing dual condensate. This feature reflects 
the conserved center symmetry: the trace of the Polyakov loop is zero and does not prefer any 
direction in the space of boundary conditions. 
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Figure 3: The integrand (L ('« + 4'V)/^ 

as a function of the boundary condition (p for two values of 
am and two temperatures, for configurations with real Polyakov loop. For non-real Polyakov loop the plot is 
shifted by ±2n/3 restoring the 27r/3-symmetry well-known from imaginary chemical potential. 
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The situation in the deconfined phase might be more confusing at first glance as the spectrum 
has a gap there: how can the vanishing chiral condensate then generate a finite order parameter Z? 
The answer is again the dependence on the boundary conditions. For boundary conditions in line 
with the original Polyakov loop, the chiral condensate persists above the critical temperature Tc. 
This has been demonstrated on the lattice first in [8] and was recently confirmed for gauge group 
SU{2) in [9] (see also [10] for a random matrix study). This mechanism ensures a finite Z and 
should actually be at work for all T > Tc. 

In the quenched case the conventional Polyakov loop is the order parameter for center symme- 
try. Under center transformations the dressed Polyakov loop behaves in the same way. Therefore, £ 
is an order parameter for center symmetry, which is underneath our numerical findings. In Ref. [1 1] 
it has been proved that all functions of D(p integrated with the Fourier factor are center symmetry 
order parameters and, what is more, have a well-defined continuum limit. 

With dynamical quarks, the center symmetry is broken by the fermion determinant. One might 
then expect a (p-dependence and a nonvanishing £ also in the low-temperature phase (and neither 
the thin nor the dressed Polyakov loop are order parameters in the strict sense). 

When evaluated in the canonical ensemble with quark number 2=1 (see, e.g., [12]^), our 
observable has an interesting interpretation. It is obtained as the derivative of the free energy in 
the 2=1 sector with respect to the mass parameter m. In the confining phase this free energy is 
infinite and thus independent of m, as suggested by our quenched analysis. Above the transition 
the free energy becomes finite and m dependent, consistent with the quenched data. 



4. Summary and outlook 

We have shown that the response of Dirac spectra to different temporal boundary conditions 
contains information about confinement. We have defined the dressed Polyakov loop £ as a novel 
deconfinement order parameter, that interpolates between the dual chiral condensate and the con- 
ventional thin Polyakov loop in the different extremes of the mass parameter m. Among other 
properties we have shown that this quantity is IR dominated and transforms non-trivially under 
center symmetry transformations (therefore, we expect the findings to be independent of the choice 
of the lattice Dirac operator). In the same spirit, many center sensitive functions of the Dirac 
operator can be defined, also in the continuum. 

In full QCD (or with the recently studied 4-fermion interactions [13]), the critical temperatures 
of deconfinement and chiral restoration could be different. It would be interesting to see what 
happens then in our formalism. As a speculative scenario, a deconfining and chirally broken phase 
would have both the dual and the antiperiodic chiral condensate finite. This is easily possible, e.g. 
if the condensate has a sine-behavior with the boundary condition. Even a phase that is confining 
and chirally symmetric could be realized, for instance by a vanishing condensate for all boundary 
conditions. 



■^We thank Philippe de Forcrand and Michael Ilgenfritz for discussions on this point. 
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